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We calculate the isospin-mixing parameter for several = —1, Tz = and — 1 nuclei from 
Mg to Sn in the particle-number conserving Higher Tamm-Dancoff approach taking into account 
the pairing correlations. In particular we investigate the role of the Coulomb interaction and the 
|T, I = 1 pairing correlations. To do so the HTDA approach is implemented with the SIII Skyrme 
effective nucleon-nucleon interaction in the mean-field channel and a delta interaction in the pairing 
channel. We conclude from this investigation that the pairing correlations bring a large contribution 
to isospin-symmetry breaking, whereas the Coulomb interaction turns out to play a less important 
role. Moreover we find that the isospin-mixing parameters for Tz = — 1 and Tz = 1 nuclei are 
comparable while they are about twice as large for Tz — nuclei (between 3% and 6%, including 
doubly magic nuclei). 

PACS numbers: 21.10.Hw21.60. Jz21.60.Cs 
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I. INTRODUCTION 

One of the most striking aspect of the structure of atomic 
nuclei is the very small violation of the isobaric invariance. 
This is so even for heavy nuclei where the Coulomb interac- 
tion is not thought a priori to act merely in ajscrturbativc 
manner. As pointed out for instance in Rcf. [11], this is due 
to the weak variation of the symmetry-breaking Coulomb 
field over the nuclear volume. It has also been suggested 
from phcnomenological and fundamental (at the levels of 
quarks having different masses) points of view that gen- 
uine isospin non-conserving parts of the strong interaction 
should be considered. They should however be rather small 
as compared to their conserving counterparts. 

As a consequence, it has been found that a nuclear 
ground state |^) may be thought as being composed of 
mostly a To = \Tz\ component where = {N — Z)/2 with 
a small Tq + 1 admixture, namely 



If) 



alTo + lT,) 



(1) 



where -I- = 1 . 

Even though in most cases the isobaric invariance may 
be flatly assumed, there are phenomena where a specific 
knowledge of the isospin mixing is needed. This is in par- 
ticular the case whenever some observed transition or reac- 
tion would be forbidden, should this invariance be exactly 
fulfilled. Interesting cases, where the isospin mixing has to 
be considered, are also related with beta-decay properties 
(see, e.g., the review of Ref. 0)- Of particular importance 
in that respect are the studies of superallowed 0"''-to-0^ 
nuclear (3 decays in the context of the tests of the CVC hy- 
pothesis (see, e.g., Ref. 3) through /t- value measurements. 
Hence, a specific determination of the effect of the isospin 
mixing is required to correct the value yielded by the crude 



isospin-multiplet approximation (determining thus the so- 
called 6c corrective term). 

Before entering, in a subsequent study, into a detailed 
assessment of the transition matrix element involved in such 
particular decays, we consider it interesting to evaluate first 
the actual importance of the isospin mixing as measured for 
instance by a^. This is the subject of the present paper. 

Presently available theoretical estimates of the isospin 
mixing fall into three different categories. 

First, one has to quote the hydrodynamical approach of 
Bohr, Damgard and Mottelson which consists in quan- 
tifying the normal modes associated with the polarization 
effect of the Coulomb field on a spherically symmetrical 
isovector density. In N = Z nuclei this approach yields the 
probability of the T = 1 component, in sole addition to 
the dominating T = component, which is given by 



= 3.5 X IQ-'^Z^A'^/^ . 



(2) 



It therefore amounts, e.g., for the Ca nucleus to about 
0.16%. In nuclei having a neutron excess, these authors 
estimate that (meaning now the probability of the \Tz \ + 
1 component over the dominating \Tz\ component) is equal 
to the value given by Eq. ^ divided by |r^| -I- 1. This 
reduction, which is expressed in terms of a factor being 
merely the square of a Clebsch-Gordan coefficient, has been 
first advocated by Lane and Soper It yields, e.g., for 
the ^*Ca nucleus, a value of of about 0.04%. 

The second class of approaches are based on shell-model 
calculations. Their success is contingent, as usual within 
such an approach, upon the relevance of the matrix ele- 
ments in use. For the description of isospin mixing, an accu- 
rate determination of Coulomb matrix elements is of course 
of paramount importance (see for instance the discussion of 
Coulomb energy differences in A = 47 and A = 49 mirror 
pairs 0). This constitutes an a priori necessary condition 
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to provide valuable answers to the question left open on the 
real importance of isospin non-conserving forces as studied 
for example to explain the isobaric multiplet yrast energies 
in Ref. [31 ■ Other concerns are related to a good description 
of radial single-particle wave functions as in, e.g., Ref. Q 
to describe asymmetry factors in parity-violating electron 
scattering. One definite difficulty of shell-model calcula- 
tions is due to the fact that they do not take into account 
any core isospin mixing, excepted of course for the no-core 
shell model calculations limited to very light nuclei (see for 
instance Ref. 0). 

One might then be inclined to think that microscopic cal- 
culations making use of phcnomcnological nucleon-nucleon 
forces should be able to describe the polarization effects of 
the Coulomb interaction at least at the mean field level, in 
a satisfactory way. Indeed, as opposed to shell-model cal- 
culations, mean-field calculations are expected to provide 
rather elaborate single-particle wave functions and they do 
not rely on any inert core approximation. However, apart 
from possible consequences of well-known symmetry vio- 
lations inherent to the mean field approximation, they re- 
quest as a next very important step to account accurately 
for the correlations. This may be done without serious a 
priori problems for RPA-type correlations, as performed for 
instance in [l^, [HI . In Ref. [l0| , it is shown that the hy- 
drodynamical ansatz of Ref. [J] underestimates the isospin 
mixing by a factor 2 to 4 (see Fig. 3 of It is im- 

portant to note that the latter approach does not include 
important correlations, namely pairing correlations. There 
are good practical reasons for such an omission. The usual 
handling of pairing correlations within a kind of Bogoli- 
ubov quasiparticle vacuum approximation as in the BCS 
or Hartree-Fock-Bogoliubov theory is totally unfit for the 
isospin mixing problem. Indeed, such an ansatz yields spu- 
rious components of both charge state particle numbers, 
giving rise in turn to a spurious mixing of T^-components 
which invalidates a priori any attempt to extract out of 
them any meaningful T-mixing properties. 

This is why we make use here of the Higher Tamm- 
Dancoff approach (HTDA) which can be interpreted as 
a highly truncated shell model built on a self-consistent 
Hartree-Fock solution [l^, [H, [13, [iBl ■ At this stage we fo- 
cus on the role of jTzl = 1 pairing correlations, which gives 
us an upper limit of the isospin mixing parameter since 
proton-neutron pairing correlations are expected to reduce 
the isospin mixing as it will be discussed below. For the 
time-being we will not evaluate the effect of RPA correla- 
tions which could be (and will soon be) easily taken into 
account into the HTDA framework. 

To determine a^, we should in principle perform a projec- 
tion of the ground state I^P) on good isospin states \TTz)- 
Assuming, however, that components higher than Tb + 1 a'l'G 
negligible, as in Eq. ([1]), we can deduce a" from the calcu- 
lation of the expectation value of the square of the isospin 
operator T in the state |\I'). Indeed, if \'^) is normalized to 
unity and assuming that the dominant contribution of the 



ground-state (GS) isospin comes from Tq — l^zl, we have 

(^'IT^I^') = (1 - a^) To (To + l) + a' (To + I) (To + 2) , 

(3) 

hence 

2 (xI/|t^|vI/)-To(To + l) 

m^) • 

The paper is organized as follows. After the derivation of 
the expression for the expectation value of in the state 
1^") in Sect. |TT1 we present in Sect. IIIII the results of the 
HTDA calculations for the GS properties and the values 
of all relevant isospin quantities, such as the expectation 
value (4'|T^|\I'), the deduced value of T and the isospin- 
mixing parameter for a large sample of nuclei. The 
main conclusions of this study are drawn in Sect. IIVI 

II. EXPECTATION VALUE OF IN THE 
HIGHER TAMM-DANCOFF APPROACH 

A. Correlated ground state in the Higher 
Tamm— Dancoff approach 

Neglecting here the proton-neutron residual interaction, 
we can write the many-body state j^*) describing the ground 
state of a nucleus as the product of the correlated states 
and 

I*) = |*("^) ® , (5) 

where, in the HTDA approach, I'l'''^-') (q = n for neutrons 
and q = p for protons) is a superposition of A^g-particlc 
Slater determinants {Ng = N for neutrons and Nq — Z for 
protons) of the form 

l*''^>=X^'^I^^^HE^i'^l*i'V (6) 

In Eq. ©, \^^^^) denotes the Hartree-Fock (HE) ground 
state and the | are n-particle-n-hole excited states 
built on |$q''^)^. The a priori complex coefficients Xo^ ^'^d 
xf'^ are determined by minimizing the energy functional 
calculated for j^I''^)). In fact, in order for the many-body 
state IvP), when constructed with real single particle wave 
functions, to be time-reversal invariant, the coefficients 
and x^rP must be real. 

B. Expression of the expectation value of T'^ 

Since is an hermitian operator (see Appendix A for 
its definition and properties), its expectation value in the 



^ For the sake of clarity in the notation, we reserve the letter 3> for a 
Slater determinant and the letter >!' for a correlated state. 
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HTDA state |*) reads 
(^IT^I*) = 



2 E xT^xfxfx^^ X 

(($(")] «(ci,W|)T2(i$(r))«i$W)) 



(7) 



Re 



where Re(z) denotes the real part of the complex number 
z. Since is a sum of one-body and two-body operators, 
the only contributions in the off-diagonal term of Eq. ([7]) are 



therefore those for which 1$ 



|1> 



, and |4"^>®|$}^'0 
differ by a particle-hole excitation of order less than or equal 

to 2. In the following it will be useful to recall that, if IfE*^"') 

and 1$''''') are Slater determinants of N and Z particles, 



respectively, then |<i>. 



1$ 



(p)\ 



is a Slater determinant of 



A = N + Z particles. Moreover, a Slater determinant 
without indication of its charge state q is to be understood 
as a product of a neutron l^j"^) and a proton Slater 
determinants. Finally we recall that the number of particles 
of each charge state q is even since we treat here even-even 
nuclei only. 

Using the expressions for the isospin operator developed 
in Appendix A, it is easy to show that the diagonal matrix 
element of can be written in the form 



-EE K^i^) 



space— spin | : 



(8) 

where the notations (fc|^) space-spin and X^^g^t") are defined 
in Appendix A. It is important to note that the sums over 
the occupied single-particle states of | $'"'') and l^^^') can- 
not be a priori reduced to sums over time-reversed partner 
states, except for the many-body states (including |<i>o)) in 
which the single-particle states are all paired. When this 
is not the case, the contributions of the form (A:|^)spacc-spin 
or space-spin (whcrc | fc) is the time-reversed partner of 
I A:)) vanish. Therefore the expectation value of in the 
Hartree-Fock ground state |<I>o) is a special case of Eq. (|5]). 



The contribution (4'|T^|5')diag of the diagonal terms in 
Eq. ([S|) to the expectation value of finally writes 



E(xl"^x?y E E \(m 



(9) 



space— spin | 



To calculate the off-diagonal matrix elements we can ex- 
ploit the fact that one of the two Slater determinants of a 
given charge state is expressed as a n-particle-n-hole ex- 
citation with respect to the other one. This gives sim- 
ple expressions for the matrix elements but for each pair 
of Slater determinants 1$,-'''') and we have to deter- 



Iji), 



mine the single-particle states 
\jn) (hole or particle states of |$o)) such that |$ 

ifiij aj^ • • -al^ttj-^ ■ ■ •Oj^l^l'''') , where ipij = ±1 is a phase 
factor determined in Appendix B. 

The non vanishing off-diagonal matrix element of in- 
volving two Slater determinants differing by a 1-particle-l- 
hole excitation ajoj with i 7^ j is given by 



~ ^kn^ln ^ ^ (-^ I ''^^) space — spin (^l^)spacc — spii 



(10) 



5kp5tp E ^^1^ 



/space— spin \"''|'^'/space— spin • 



\k)s 



Finally, for two Slater determinants differing by a 
2-particle-2-hole excitation ^ '^ji with {11,12} H 
{ji,j2} = 0, we have 

{^km Sk2P ^lip ^£2" + ^fclP ^k2n ^Im ^e.2p) X 

(£1 |fci)spacc— spin (^2 | ^2) space— spin (H) 
— {5kip Sk2n ^£ip Si^n + Skin Sk2p ^^m ^fap)^ 
(^1 |^2)spacc— spin (^2 | ^l) space- spin • 

From Eqs. PH)) and pT|) we deduce that the non van- 
ishing off-diagonal contribution (4'|T^|4')off-diag to the ex- 
pectation value of takes the form 



i.j i' = lplli(i) 

+ 2EU!"Vxf E x?(($r^i«($f i)T^(i$!"')«i$?)) 

i,j i' = lplli(i) 

/=lplh(j) 
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The first two terms of (^'|T^ |5')off-diag are calculated using 
respectively the first or the second term of the right hand 
side of Eq. ([10]), whereas the third term of (^'|T2|5')off-diag 
corresponds to one of the four series of S products in the 
right hand side of Eq. ([TT|). In practice, the sum of all off- 
diagonal terms is at least two orders of magnitude smaller 
than (*|t2|*)diag. 



C. Limiting cases 

We consider in this subsection two relevant limiting cases: 
the Hartree-Fock limit and the limit of identical neutron 
and proton single-particle states. We will refer to theses 
limits in Sect. Ill to interpret some results. 

In the Hartree-Fock limit where x^-^' = SiQ, we can de- 
duce from Eqs. ([S]), ^TU\\ and (fTT|) that the off-diagonal ma- 
trix elements vanish. The expectation value of thus 
simply becomes 



A {N- Zf 



2+ 4 



space— spin I 



In the limit where the neutron and proton single-particle 
states are assumed to be identical, the diagonal contribu- 
tion (O to (*|T2|^^ becomes 



(vI/|T^|vI/)di,g =ro(To + l) 



(ij)#(0,0) 



(14) 

where Tq = \Ty\ and 7^($|"^ $^-^^) denotes the relative ex- 
citation order of | with respect of (see Eq. (|B-8|I 
of Appendix B). Therefore the isospin-mixing parameter 
takes, in this model case of identical neutron and proton 
single-particle states, the simple form 



1 



2(ro + 1) 



E (xfVfT7^($("^ (15) 



(ij)#(0,0) 



III. RESULTS AND DISCUSSION 

We study the isospin symmetry breaking through the 
isospin-mixing parameter o? defined in Eq. ([4]) for = —1, 
= and = 1 nuclei of eight elements, namely Z = 12 
(Mg), Z = 16 (S), Z = 20 (Ca), Z = 24 (Cr), Z = 28 (Ni), 
Z = 36 (Kr), Z = 40 (Zr) and Z = 50 (Sn). 

To evaluate the expectation value (4'|T^|5') we need a re- 
liable description of the ground states of these nuclei. For 
that purpose, we follow the two-step approach of Ref. [lH| 
where it was applied to study GS pairing properties of 
N = Z nuclei in the mass A ^ 70 region. Since the above 
considered nuclei exhibit no triaxial deformation in their 



ground state, we can search for GS solutions possessing 
axial symmetry. In the first step, we determine the GS 
deformation within the Hartree-Fock-BCS (HFBCS) ap- 
proach. To do so. we use the Skyrme interaction in its SIII 
parametrization [l6| in the mean-field channel, and the se- 
niority force in the pairing channel. For the latter we re- 
tain the same set of parameters as in Ref. [l^ , where they 
were adjusted to reproduce experimental odd-even mass dif- 
ferences through a 3-point formula (Gq"^ = 17.70 MeV, 

G[,^' = 15.93 MeV, Ae 6 MeV and /i = 0.2 MeV). In 
practice, we use 15 oscillator major shells to expand the 
single-particle states on the cylindrical harmonic-oscillator 
basis and optimize the basis parameters at the GS deforma- 
tion so as to obtain the lowest HFBCS binding energy. In 
the second step, we calculate GS properties in the HTDA 
approach from the above HFBCS solution. The residual 
interaction employed is the delta interaction of Ref. [l^ 
adjusted in the same way as above for the seniority force 
but with Ae = 12 MeV. The optimal values V^''^ of the 



strength were found to be Vq 



(") 



-340 MeV.fm^ and 



Vr 



(p) 



—306 MeV.fm'^ (this fit has been performed on 
(13) the neutron pairing strength upon the simple approxima- 



tion that Vq^'^ is quenched by 10% with respect to Vq 
because of the anti-pairing contribution of the Coulomb in- 
teraction). However, keeping the same interaction strength 
throughout the whole considered nuclear region, we have 
taken care of the well-known A^^^^ energy scale by varying 
the active pairing window: Ae = 12 x (72/A)^/^ MeV and 
H = 0.2x (72/A)i/3 MeV, which yields for A = 72 the same 
window parameters as those of Ref. [l^ . 

The GS properties calculated here are the charge radius 
Tc, the (32 deformation parameter (see Appendix C), the 
mass quadrupole {Q20) and hexadecapole (Q40) moments. 



(") 



the neutron and proton pair-condensation energies 



(9) 



- p) , which is equal to the 
^ and Vi = ^/pii, with p 
Tsl). and the total 



the trace of the operator p{l 

sum J2i '^i'^i with Ui = \/l — v:^ 
being the one-body density (see Ref. 
binding energy E^,. The results are reported in Table [H 

The resulting HTDA ground state j^*) is then used to cal- 
culate the expectation value of the operator. In practice 
the off-diagonal term ([T^ turns out to be negligible with 
respect to the diagonal contribution ([9]) and therefore can 
be safely omitted in the calculations. Then, from the value 
of (*|T2|^'), we deduce the T- value defined by 



(^IT^I^-) = T{T+l) 



(16) 



We present the values of (^|T^|5'), T and in the columns 
labeled "HTDA" in Table M and in Fig. [1] we show the 
variation with Z of the isospin-mixing parameter within 
the HTDA approach for the above twenty four nuclei. 
Apart from a dip around ^^Ni, increases with Z, faster 
for the N = Z nuclei than for the others. Moreover the 
dip is deeper for the former nuclei. We also note that the 
isospin-mixing parameters for = — 1 and ~ I nuclei 
are very similar and are about a factor of two smaller than 
for Tz = nuclei. 

A precise determination of the isospin-mixing parameter 
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TABLE I: Ground-state properties of the twenty-one studied nuclei calculated within the HTDA approach, from left to right: 
the charge radius r^, the deformation parameter (calculated as in Eq. (|C-6|) of Appendix C), the mass quadrupole {Q20) and 
hexadecapole (Q40) moments, the neutron and proton pair-condensation energies -E'^^^j. 



Nucleus r, (fm) Q20 (fm') Q40 (fm*) ^^''^^ (MeV) 
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requires that some great care be exerted in the calculations. 
This is illustrated in some typical examples in Appendix D. 
Here we merely discuss the most important points. 

First of all, we need to make sure that we have obtained 
a perfect consistency between the wave functions and the 
mean field including its Coulomb isospin-brcaking part. In 
Appendix D wc show that a poor convergence of the iter- 
ative process may lead to drastic distortions in the isospin 
mixing evaluations. 

A second important point is related to the quantal char- 
acter of the assessed quantity. As exemplified in Ap- 
pendix D, it appears that in order to get reliable val- 
ues, one should include almost all Slater determinant com- 
ponents 1$,;) of the correlated wave function even 
those appearing in j^*) with a relatively minute probability 

{Xi^^)'^i because of constructive interference effects. 

A last technical point is worth noting here. It deals with 
the question of the independence of our results with the 
harmonic oscillator basis parameters b and q (see their def- 
inition, e.g., in Ref. [l3|) in the expansion of the single- 



particle wave functions. One might have been concerned by 
the fact that this optimization has been performed merely 
at the level of the preliminary HFBCS calculations and not 
at the final stage of our HTDA approach. However it has 
been checked that an energy optimization of HTDA results 
leaves unchanged the calculated values, as shown on one 
example in Appendix D. 

We now investigate several sources of isospin symmetry 
breaking. 

A. Roles of the neutron-proton mass difference and 
the Coulomb interaction 

We investigate separately the sensitivity of our results to 
the neutron-proton mass difference and to the presence of 
the Coulomb interaction. 

It turns out that the former plays virtually no role at 
all. For ""^Ca, for instance, upon suppressing the one-body 
center of mass correction (involving a 1/A term, ambigu- 
ous in this context), we found that the mass difference is 
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TABLE II: Expectation value of T^, isospin T from Eq. (|16|) and isospin-mixing parameter from Eq. Q calculated within 
the HTDA and HFBCS approaches at the ground states determined in Table |l] The columns labeled "HP" correspond to the 
contributions to the above three quantities coming from the Slater determinant |<J>o) in the HTDA ground state expansion of 
Eqs. ([5l and (|6]). The values given in italic are obtained without Coulomb interaction. 
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FIG. 1: Variation with Z of the isospin-mixing parameter 
calculated within the HTDA approach for Tz = —\ (open cir- 
cles), T, = (open triangles) and T~ — 1 (open squares) nuclei. 



responsible for a variation in of the order of one part in 
10^. It is neglected in the remainder of the paper. 

As can be seen from the columns "HTDA" of Table 
where the results obtained without the Coulomb interaction 
are shown in italic, the Coulomb interaction has a more 
noticeable effect which still remains rather weak for the 
light nuclei considered here. In this comparative study, we 
do not include the results for "'^Cr, ^"'Kr and ^^Zr because 
the GS solutions without Coulomb interaction differ too 
much from the ones obtained in the full calculations to make 
a comparison meaningful. 



B. Role of the particle- number conservation 

A very interesting issue consists in investigating the qual- 
ity of the particle-number conserving pairing treatment 
(here in the \Tz\ = 1 channel only) obtained by the HTDA 
approach as compared to approximations (as in the HFBCS 
calculations preliminary to our HTDA evaluation) which 
violate particle-number conservation. For that purpose, 
we evaluate the expectation value of from the HFBCS 
ground state and deduce the value of the isospin-mixing 
parameter through Eq. (|4]). The expectation value of in 
a BCS state normalized to unity, noted |BCS), reads 



1 



(BCSIT^IBCS) = A+^(N 



2 ^"''"'^f '"^'^'^k I (* I fc) space-spin 



i>0 
|2 



(17) 



j>0 



k>0 



states of the form respectively. The resulting val- 

ues of (BCS|T^|BCS), T and are reported in the columns 
labeled "BCS" of Table [H. 

In cases where pairing correlations are ineffective in the 
BCS treatment (one is then below the phase transition to 
the superfluid phase), the value of is very small, since in 
that case essentially a single Slater determinant is describ- 
ing the nuclear state and the particle number is trivially 
conserved. In contrast, for those nuclei where pairing plays 
a non negligible role, the values of the parameter turn 
out to take on completely unrealistic values as, e.g., for the 
nuclei "^^Ca or ^''Zr. In fact, as shown in Fig. [21 there is 
a strong correlation between calculated in the HFBCS 
approach and the particle-number fluctuation AA^ -|- AZ in 
the BCS state, where AA^^ is defined by 



AN, = J(BCS|Af2|BCS) 



Ar2 
q 



(18) 



With the exception of ^^Kr and ^"Zr, all the points lie ap- 
proximately on a straight line in the (a^ , A A'' -|- AZ) plane 
as can be seen in Fig. [2l 




where the sums J2i>o '"^ Sfc>o ''^^ ^^'^ ^i>o '^^'■^^ 
tron, proton and all pairs of time-reversed single-particle 



0.5 1 1.5 2 2.5 3 3.5 4 
AN+AZ 

FIG. 2: Correlation between and the particle-number fluctu- 
ation AN + AZ within the HFBCS approach for all twenty four 
Tz = -1, Tz =0 and Tz = 1 nuclei. 



C. Role of the pairing correlations 

In order to assess the importance of pairing correlations 
on the isospin-mixing rate, we also calculate the expec- 
tation value of in the Slater determinant |$o) using 
Eq. The resulting values for ($o|T2|$o>, T and 

are reported in the columns labeled "HE" of Table [TTl In 
general ($o|T^|$o) is different from the value that would 
result from a pure HE calculation because |$o) is the Slater 
determinant built up from the single-particle states result- 
ing from the first-step HFBCS calculation. This difference 
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vanishes of course for nuclei in which BCS predicts no pair- 
ing correlations, which is the case here for the doubly magic 
nuclei as well as ^"'Mg and '^^S. 

For nuclei exhibiting weak pairing correlations, the "HF" 
results are, quite expectedly, close to the HTDA predic- 
tions. Otherwise, the HTDA results are significantly larger 
than the "HF" ones. The = pairing correlations arc 
therefore an important source of isospin symmetry break- 
ing. This is conspicuous from Fig. [3] which shows the strong 
correlation between the variations of and UiVi with 
Z for = Z nuclei. 



is simply given, here, by 



6 

._ 5 
>_ 




10 15 20 25 30 35 40 45 50 55 60 
Proton number Z 

FIG. 3: Correlated variations with Z of (solid circles) and 
UiVi (open circles) for = nuclei calculated within the 
HTDA approach. 



(") ^ip)\ _ 



7^($^^^^) 



if i = 0, j = 0; 

2 if i = 0, j ^0 or i^O, j = 0; (20) 
4 otherwise. 



Since the states j^E''^'^^) are normalized to unity, we have 

(21) 



and we finally obtain 
1 



Tn + l 



2-(x^"T-(xr)' 



(22) 



Since Eq. (|22p overestimates the importance of one-pair ex- 
citations in 1^) through the estimates of Eq. (|20[1 and given 
the very small contribution to j^*) coming from the particle- 
hole excitations other than one-pair excitations, we con- 
clude that the value of calculated with Eq. (QUI should 
lie between the values obtained in the full HTDA calcula- 
tions without and with Coulomb interaction. In the case 
of "^OCa, we find = 0.9804 and {xo^f = 0.9825. 

This yields « 3.7%, which is slightly larger than the 
3.5% obtained in the HTDA calculation without Coulomb 
interaction and smaller than the value of 4.2% from the full 
HTDA calculation, as expected. 

Finally, it is interesting to note (see Table |Tl| that the 
value of for a given nucleus obtained in a full HTDA 
calculation can be written, to a good approximation, as the 
sum of the "no Coulomb" HTDA result (including pairing 
correlations) and the "HF" result (no pairing, but including 
the full Coulomb field). 



The large values of the isospin-mixing parameter found 
particularly in the = doubly-magic nuclei can be ex- 
plained as follows taking the example of *°Ca. As can be 
seen in Table [Hi the Coulomb contribution to the isospin 
mixing in HTDA calculations can be considered as small. 
We can therefore make the approximation that the neu- 
tron and proton single-particle states of '"'Ca are virtually 
identical and use Eq. p3|) to estimate a^. Moreover, the 
largely dominant contributions to the particle-hole excita- 
tion expansion of j'l') in Eqs. ([5]) and ([6]) come from one-pair 
excitations so that we can write approximately as 



1 



2(To + 1 



(19) 



where the relative excitation order of with respect to 



D. Discussion 

In Fig. |4] we compare, for the above eight N — Z nu- 
clei, the values calculated in our HTDA model with the 
estimates obtained in the hydrodynamical model of Bohr, 
Damgard and Mottclson M, and with the calculations by 
Hamamoto and Sagawa in the Hartree-Fock-plus-RPA 
approach with the SHI Skyrme interaction. Each model 
predicts an increasing trend of with Z. The HTDA 
approach, as presently applied with \Tz\ = 1 pairing corre- 
lations only, predicts a larger isospin mixing than the RPA 
calculations (which do not include pairing correlations). 

It is important to recall that only jT^j = 1 pairing corre- 
lations arc considered here. We expect the values of ob- 
tained by including, in addition, pairing correlations in the 
= channel to be smaller than the present values. In- 
deed, in presence of proton-neutron correlations, the \Tz \ = 
1 pairing correlations, being effected by a smaller proba- 
bility amplitude, would contribute less to the total HTDA 
wave function. Correlatively, they would be replaced essen- 
tially by configurations of the type aJaj|<I>Q"^) (E) a\a(\^'^^) 
where the neutron and proton hole states \j) and \£) on 
the one hand, the neutron and proton particle states \i) 
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FIG. 4: Variation with Z of the isospin-mixing parameter 
of T2 = nuclei calculated within the HTDA approach (solid 
circles), by Hamamoto et al. [l^l in the RPA approach (open 
circles) and by Bohr et al. 0] in a hydrodynamical model (dashed 
line) . 



and |fc) on the other hand, are similar. The relative excita- 
tion order between two such one-particle-one-hole neutron, 
one-particle-one-hole proton configurations would thus be 
on average smaller than that between two one-pair exci- 
tation neutron, one-pair excitation proton configurations 
which reaches about 4 from Eq. (|20p . However, the further 
addition of RPA correlations, which can be treated in the 
HTDA framework on the same footing as the pairing corre- 
lations, would compensate the effect of the = pairing 
correlations and the net result may be close to the present 
result. 



IV. CONCLUSION 

We have shown that the HTDA approach is a reliable 
model to address the isospin-mixing issue because such an 
approach can include the pairing correlations in a consis- 
tent way while conserving the particle-number, in contrast 
to Hartree-Fock-Bogoliubov and Hartree-Fock-BCS treat- 
ments. 

From HTDA calculations, we have learned that the dif- 
ference between the neutron and proton masses has a negli- 
gible impact on the expectation value of and that the ef- 
fect of the Coulomb interaction is rather small as compared 
to the effect of the \Tz\ = 1 correlations. Moreover the 
stronger isospin symmetry breaking is found in the N = Z 
nuclei. 

To obtain a more complete description, both neutron- 
proton pairing and RPA-type correlations need to be taken 
into account. Both of these can be included in the HTDA 
framework in a consistent way. It is expected that these two 



types of correlations affect the present results with opposite 
signs. Such a study is currently under way. 
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APPENDIX A. ISOSPIN OPERATOR AND 
ONE- AND TWO-BODY MATRIX ELEMENTS 

The operator can be written as the sum of a one-body 
operator Oi and a two-body operator O2 acting in the Fock 
space 



O2 



i 

= i^2(t,® 



(A-1) 
(A-2) 

(A-3) 



where t,- 



t'r. 



ix -\- ty ® iy + iz ® iz- We define the 
one-body and two-body operators di and 62 acting in the 
one-particle space and the two-particle space, respectively, 

by 



01 = t^ 

62 = 2 (t (g) t) = 2 (4 (g) 4 + iy (g) iy + iz (g 4) 

Introducing the operators and i_ defined by 

i+ =ix + ii'u 



1-1/ , 



t. 



it: 



V ' 



we can rewrite 62 as 

62 = i+ (g t_ + £_ (g £+ + 2 (4 (g 4) 



(A-4) 
(A-5) 



(A-6) 
(A-7) 

(A-8) 



The Hartree-Fock basis is built up from the single- 
particle states generically noted They describe either a 
neutron state or a proton state, so they are eigenstates of 
the isospin operators and tz 



tz\i) = n\i) 



neutron 



proton . 



(A-9) 
(A-10) 



In practice, we expand the single-particle states \i) on the 
cylindrical harmonic oscillator (HO) basis {|q;)} as follows 



(A-ll) 
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where a stands for the 4 quantum numbers riz, nj_, A 
(eigenvalue of £z) and S (eigenvalue of s^), t — 1/2 and 
Tj = ±1/2 depending on the nature of the particle. 

The action of the operators t_ and t+ on the single- 
particle states \i) is given by 



where \ij) 



In) 



(A-12) 
(A-13) 



The matrix elements of t_ and t+ in the Hartree-Fock basis 
thus write 

{i\i^\k) = (5ip4n(i|fc>spaco-spin , (A-14) 
{i\i+\k) = (5m 4p(i I ^) space-spin • (A-15) 

In the above equations and elsewhere in this paper, the sub- 
script "space-spin" attached to an overlap of single-particle 
states means that the overlap is restricted to the space 
and spin variables only. This allows to consider such over- 
laps between two nucleonic states corresponding to different 
charges. Since the time reversal operator does not act on 
isospin, we have 

{i\i±\k)=0, (A-16) 
where |A:) is the time-reversed conjugate state of \k), and 
(i\i±\k) ^ {i\i±\k) . (A-17) 
From Eq. (jA-9|) we easily get 



(A-18) 



and using Eqs. (|A-8p . (|A-14p and (jA-15|) . we can write the 
two-body matrix element {ij\62\k£) as 

(^l^')spacc— spin (j |^) space— spin ■ 

To close this appendix, we recall useful expressions for 
the expectation value of one-body and two-body operators 
in a Slater determinant and related matrix elements. 
For a one-body operator Oi we have 

(<i>,|Oi|<I>,) = ^(fc|6i|fc), (A-20) 

md,aia,\^,)^S^Sim\k), (A-21) 
where the sum ^ runs over the occupied single-particle 

fcG*i 

states |fc) of the Slater determinant In Eq. (jA-21[) 

and below, Sj (resp. S^) is equal to 1 if |j) (resp. \k)) is a 
hole state (resp. particle state) with respect to |3>i) and 
otherwise. For two-body operators we have 
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jk\02\jk) , 



(A-22) 
(A-23) 



($.|024,«l«.l«.J*.> - {Jlj2\d2\klk2) , 

(A-24) 



APPENDIX B. COMPARISON OF TWO SLATER 
DETERMINANTS 

Let us consider two Slater determinants and 
built from the same set of orthonormal single-particle basis 
states. They may therefore be thought of as n-particle- 
n-hole and n'-particle-n'-hole excitations on a reference 
Slater determinant which may be chosen as the Hartree- 
Fock ground state |<I>o) 

4i ■ ■ ■ "L^f-i ■ ■ ■ "fc" i^o) , (B-i) 

|$j)==at^---at^,aci---ac„, |$o), (B-2) 

where /3i < • • • < /3„ and 71 < • • • < 7„' are two sets of 
particle states (with respect to |$o)) and bi < ■ ■ ■ < bn and 
ci < • • • < c„' are two sets of hole states. We can express 
|(f>j) as a function of as 

= flt^ • • - ajj, ,aci • ■■ac^,al ■ ■ ■ al^ap^ ' ' ' a/3i l^i) ■ 

(B-3) 

We denote by Ti the set of hole states in common between 
{bi, ■ ■■ ,6„} and {ci, • • • ,c„'} 



n= {ch,=bh^,li^ is:Nh} , 



(B-4) 



where A^^ is the number of hole states in common. Sim- 
ilarly V is the set of particle states in common between 
{/3i, • • • ,/3„} and {71, • • • ,7„'} 



^={7p; =l3p,,l^i^Np} 



(B-5) 



where Np is the number of hole states in common. There- 
fore it can be shown that 



V. n< n 



k=l 



k—7i 



k=l 



n ^'^O'* 



k—n 



where the associated relative phase is given by 



= (-1) 



(B-6) 



(B-7) 



Changing the order of the creation and/or annihilation op- 
erators in Eq. (|B-6p would change the sign of ipij . 

Finally the relative excitation order 7?.($i, <&_,) between 
the two Slater determinants and defined as the 

number of creation (or annihilation) operators in Eq. (|B-6p , 
is simply given by 



7^($,, $j) = n + n' - [Nh + iVp) 



(B- 
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APPENDIX C. NUCLEAR SHAPE AND SIZE 
QUANTITIES 

Starting from the nuclear shape determined in a self- 
consistent way by the HTDA solution, we can extract 
a quadrupole deformation parameter (32 by approximat- 
ing the nuclear shape by the equivalent spheroid hav- 
ing the same root-mean-square mass radius and mass 
quadrupole moment Q20 as the actual nucleus. The semi- 
axes c (along the symmetry axis) and a (in the perpendic- 
ular direction) are related to and Q20 through 



with 



Arl^= J <fvp{v)v^ = i(2a- 
Q2o = 2 j d\p{r)r^P2{cos0) = ^ ^ (c' 



2 ■ c^) 



(C-1) 
') , (C-2) 



where A = N + Z, p(r) is the isoscalar nuclear density 
(sum of neutron and proton contributions) and P2 is the 
Legendre polynomial of degree 2. 

The [32 parameter is then calculated for this equivalent 
spheroid by expanding the nuclear radius in polar coordi- 
nates according to the /3i-parametrization [l8l | 



m 



V 1 — a cos^ 9 

/ 00 N 



with 



1=1 



1-^ 

C2 



(C-3) 
(C-4) 

(C-5) 



This allows us to derive the analytical expression of /32 for 
the equivalent spheroid as a function of a as 



_3_ / -j^ _ 1 _ ^ 

2a \ Arcsin ^/a 



/5tt 



( 1 



In (\/ — a+Vl^ct) 



a e]0;l[ 
a = Q 



(C-6) 

As for the mass hexadecapole moment (540, we calculate 
it using the following expression with usual notation 



Ho 



d'vp{v)r^Y^{e) 



(C-7) 



Finally, the charge radius Tc is calculated as in Refs. [IE 
[131 through 



(C-8) 



where Pp(r) is the proton density and /p(x) denotes the 
proton form factor. With a Gaussian form for the latter, 
/p(x) exp{-x'^/r^)/{ro ^/7^), we have 

rl=rl + lrS, (C-9) 



(C-10) 



In our calculations we choose to use the value tq = 0.65 fm 



(|r 2 = 0.64 fm^) from Ref. [2C 



APPENDIX D. 



TECHNICAL ASPECTS OF THE 
CALCULATIONS 



We illustrate in this appendix the importance of several 
technical aspects of the mean-field calculations that can 
have a substantial impact on the quality of the results for 
the isospin-mixing parameter . 

First, the choice of the reference Slater determinant |$o) 
for HTDA calculations is very important. Figure [5] shows 
the variation of the isospin-mixing parameter calculated 
within the HTDA approach with the number of prelimi- 
nary HFBCS iterations, in the case of ^^Cr. As an initial 
potential we choose the Woods-Saxon potential including 
a spin-orbit term with the same parameters for neutron 
and protons (hence without Coulomb interaction for the 
first iteration). From the decreasing and saturating trend 




20 40 60 80 100 
Number of preliminary HFBCS iterations 



FIG. 5: Variation of the isospin-mixing parameter calculated 
within the HTDA approach with the number of preliminary HF- 
BCS iterations, in the case of **Cr. 

obtained here, we conclude that the consistency of the un- 
derlying mean-field from which |$o) is determined plays an 
important role. In other words a poor mean-field is a source 
of spurious isospin symmetry breaking. 

The second aspect of importance is the convergence of 
with the number m of Slater determinants contributing 
to 1*) retained in the calculation of (^iT^j*). The Slater 
determinants li*-^^) entering the expansion ^ of are 



12 



The results are shown 



arranged in decreasing order of \xl 
in Fig. [HI whereas in Fig. [7] we present the variation of 
the neutron and proton Slater determinant amplitudes \xm\ 
with TO. 



< 




500 1000 1500 2000 

Number m of components |<I>|> 



Finally, we have checked for three nuclei that the values 
of the isospin-mixing parameter calculated in the HTDA 
approach is not sensitive to the procedure of optimization 
of the harmonic oscillator basis parameters b and q (with 
the notation of Ref. namely at the level of the prelimi- 
nary HFBCS calculations or at the final stage of our HTDA 
approach. The optimized values of h and q obtained in each 
of the two schemes are reported in Table IIIII together with 
the resulting values of for ^''Mg, ^°Ca and ^^Zt. We 
recall that the cylindrical harmonic oscillator basis used in 
all calculations contains 15 major shells, which corresponds 
to iVo = 14 in the notation of Ref. pTf . 



FIG. 6: Variation of with the number of Slater determinants 
from the expansion of the HTDA ground state \^), that 
are used to calculate . 




500 1000 1500 2000 



Number m of components |<I>|> 



TABLE III: Optimized values of h and q obtained in the two 
optimization schemes (at the level of the preliminary HFBCS 
calculations or at the final stage of the HTDA calculations) and 
resulting values of the isospin-mixing parameter calculated 
in the HTDA approach for ^''Mg, ""Ca and **°Zr. 



Nucleus 


HFBCS optimization 
h q o? (%) 


HTDA optimization 
h q a" (%) 


24Mg 


0.65 1.28 3.001 


0.65 1.05 2.999 


«Ca 


0.66 1.00 4.228 


0.66 1.00 4.228 


»«Zr 


0.60 1.37 5.229 


0.58 1.45 5.264 



FIG. 7: Variation of the neutron (open circles) and proton (open 
triangles) Slater determinant amplitudes \xm\ with the number 
m of Slater determinants retained from the expansion of l^*) for 
the calculation of o? . The first point is for m — 1, the dominant 
component in \^). 
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